The universal behavior of one-dimensional, multi-species branching and annihilating 

random walks with exclusion 
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A directed percolation process with two symmetric particle species exhibiting exclusion in one 
dimension is investigated numerically. It is shown that if the species are coupled by branching 
(A — * AB, B — > BA) a continuous phase transition will appear at zero branching rate limit belonging 
to the same universality class as that of the dynamical two-offspring (2-BARW2) model. This class 
persists even if the branching is biased towards one of the species. If the two systems are not coupled 
by branching but hard-core interaction is allowed only the transition will occur at finite branching 
rate belonging to the usual 1 + 1 dimensional directed percolation class. 
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The study of phase transitions in low dimensions is 
an interesting and widely investigated topic [jl]||] (since 
the mean- field solution is not valid). The research of 
non-equilibrium phase transitions occurring in one di- 
mensional coupled systems has drawn interest nowadays 
|p[-|T^|. Several models have been found with transitions 
that do not belong the robust directed percolation (DP) 
class ]l5|-p7| or to the parity conserving (PC) class [ ^8|Jl9| 
which are the most prominent ones among one compo- 
nent systems. Particle blocking which is common in one 
dimension has not been taken into account in field theo- 
retical description of these models yet [^0|,|l] . It has been 
known for some time that the pair contact process ]24| ] 
can be regarded as a coupled system that exhibits DP 
class static exponents while the spreading ones depend 
on initial densities |p5|| . The field theoretical investiga- 
tion of Janssen predicts that in coupled DP systems 
the symmetry between species is unstable and generally 
a phase transition belongs to the class of unidirectionally 
coupled DP where coupling between pairs of species is 
relevant in one direction only. Such systems have been 
shown to describe also certain surface roughening pro- 
cesses H,|n|. 

Recently we have shown |22[ that in the two- 
component annihilating random walk (A A — > 0, BB — > 
0) owing to the hard-core interaction of particles dynam- 
ical exponents are non-universal. Some consequences of 
hard core effects for random walks in one dimension have 
been known for some time already |P3|] . 

Very recently simulations gave numerical evi- 

dence that in the two-component branching and annihi- 
lating random walk (2-BARW2) the lack of particle ex- 
change between different species results in new universal- 
ity classes in contrast to widespread beliefs that bosonic 
field theory can well describe these systems. The critical 
exponents obtained numerically suggest that the location 
of offspring particles at branching is the relevant factor 
that determines the critical behavior. In particular if the 
parent separates the offsprings: 



1) A — * BAB the steady state density will be higher 
than in the case when they are created on the same site: 

2) A — ► ABB for a given branching rate because in 
the former case they are unable to annihilate with each 
other. This results in different order parameter expo- 
nents for the symmetric (2-BARW2s) and the asymmet- 
ric (2-BARW2a) cases (J3 S = 1/2 vs (3 a = 2 for 1) and 2), 
respectively) . 

Hard-core effects are conjectured to cause a series of 
new universality classes in one dimension jfj]. In this 
paper I point out that probably only a few universal- 
ity classes emerge as the consequence of particle exclu- 
sion, because other symmetries and conservation laws 
(like that of the PC class) will become irrelevant. 

In the present study first I show that in case of the two- 
component single off-spring BARW model (2-BARW1) 
defined as 
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a continuous phase transition will occur at zero branch- 
ing rate limit (a = 0) like in the 2-BARW2 model where 
they are equivalent and therefore the exponents on the 
critical point must be the same as those determined in 
p,R0|. Furthermore I show that the order parameter 
exponent describing the singular behavior of the steady 
state density near the critical point coincides with that 
of the 2-BARW2s model. 

The particle system was simulated on a lattice with 
size L = 4 x 10 and periodic boundary conditions for 
different er-s (with A = d = 1 — a condition). The ini- 
tial condition was uniformly random distribution of As 
and Bs with a total concentration 0.5. The evolution 
of the density was followed until steady state has been 
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reached plus t ~ 10 4 Monte Carlo sweeps (throughout 
the whole paper t is measured in units of Monte Carlo 
sweeps (MCS) of the lattices). As Figure El shows a phase 
transition occurs at a a — &b = indeed. 
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FIG. 1. Steady state density as the function of a in 
the one-dimensional 2-BARW1 model. Circles correspond 
to pa + Pb for oa = &b\ crosses to pA, stars to pB when 
&A = tb/2. 

The order parameter exponent has been determined 
with the local slope analysis of the data 
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providing an estimate for the true asymptotic behavior 
of the order parameter 



/3 = lira p eff (a). 



(7) 



As one see on Figure || /3 e // extrapolates to P = 0.50(1) 
with a strong correction to scaling like in case of the 2- 
BARW2s model |?J . The coincidence of this off-critical 
exponent in addition to the equivalence of processes at 
the critical point assures that they belong to the same 
universality class. 
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If we destroy the symmetry between species by the 
branching rates: a a = cs/2 we still get the same or- 
der parameter exponents (P — 0.50(1)) for both species 
(Fig.^). Therefore this universality class is stable with 
respect to coupling strengths unlike the coloured and 
flavoured directed percolation 

It is also insensitive whether or not the parity of parti- 
cles is conserved meaning that the A — > BAB process can 
be decomposed to a sequence of A — > AB, AB — > BAB 
processes. This may seem to be quite obvious when par- 
ticle exchange is not allowed, and if locality is assumed. 
By the choice of parameters d = 1 — a in the neighbour- 
hood of the critical point the the diffusion is strong and 
the locality condition is not met. Still the two process 
share the same critical behavior. 

If we decouple the two systems and allow hard-core 
exclusion only: 



(8) 
(9) 

(10) 

(11) 
(12) 



the critical point will be shifted to a — 0.81107(1) and 
DP like density decay can be observed on the local slopes 
defined as 
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(where we use m = 8 usually) (see Figure 
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FIG. 2. Effective j3 in the 2-BARW1 model as the function 
of <r ' 5 . Different symbols denote the same as in Fig.jj]. 



FIG. 3. Effective decay exponent a e ff as the function of 
1/t in the decoupled two component DP model. The system 
size is L = 4 x 10 4 , the decay is followed for 2 x 10 5 MCS. The 
different curves correspond to: a = 0.81103, 0.81107, 0.81109, 
0.8111 (from bottom to top). 
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One can not observe any relevant correction to scaling 
here, the most straight curve corresponding to the criti- 
cal one (cr = 0.81107) extrapolates to a = 0.158(2) which 
agrees very well with the f3/v\\ = 0.159464(6) value of the 
1 + 1 DP class value that can be found in the literature 
p6[ . This is different from the case of coupled annihi- 
lating random walk, where the blocking causes marginal 
perturbation to the standard decay process [^2| . 

One can generalize the results by taking into account 
that neighbouring AA and BB offsprings decay very 
quickly and therefore irrelevant for the leading scaling 
behavior. 

Conjecture : In coupled, one dimensional N- component 
BARW systems with particle exclusion and branching 
processes like: A -> BABB, A -> BAAA, A -> BAC ... 
leaving behind non-reacting neighbouring particles which 
block each other the universality class of a phase tran- 
sition will be the same as that of l-BARW2s. If the 
branching creates only pairs that can annihilate imme- 
diately (like: A — > BAAB . . . etc.) the class of transition 
will be the same as that of the 2-BARW2a model. We can 
also conclude that in case of reaction-diffusion processes 
where spontaneous decay is allowed: 2A — > A, A — > the 
blocking effect between dissimilar species is irrelevant. 

It is very likely that the transition of a very recently 
introduced ladder model || also belongs to this class. 
This model is composed of two one dimensional subsys- 
tems following BARW at the critical point and coupled 
by ladder links. In the supercritical region by updating 
an active site one can create an offspring on the other sub- 
system or increase the inactivity level of that site. For 
small coupling strength (s = 1) the very few blocking 
events can not introduce relevant blocking on the other 
subsystem and the scaling exponents agree with those of 
the coupled BARWe model without exclusion (20). For 
stronger coupling strength (s — 2) there are more block- 
ing possibility resulting in l-BARW2s scaling exponents. 

In conclusion I have shown that the one dimensional 
two species coupled BARW with exclusion and one off- 
spring has the same critical transition point as that of the 
2-BARW2s model investigated earlier. The hard-core in- 
teraction itself is not sufficient to cause deviation in scal- 
ing behaviour from that of DP. A conjecture is given with 
regard the universality classes in coupled BARW systems 
exhibiting particle exclusion. 
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